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We show that non-reciprocal bands can be formed in a magnetized periodic chain of spherical 
plasmonic particles with two particles per unit cell. Simplified form of symmetry operators in dipole 
approximations are used to demonstrate explicitly the relation between spectral non-reciprocity and 
broken spatial-temporal symmetries. Due to hybridization among plasmon modes and free photon 
modes, strong spectral non-reciprocity appears in region slightly below the lightline, where highly 
directed guiding of energy can be supported. The results may provide a clear guidance on the design 
of one-way waveguides. 
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I. INTRODUCTION 

Breaking Lorentz reciprocity mm in optics has been 
of great interest to physicists for many decades. In re¬ 
cent years, the asymmetry in dispersion relation k) ^ 
cj(k) (i.e, spectral non-reciprocity) [3, [4] has drawn a 
lot of interest because of its possible topological na¬ 
ture [5 and potential applications such as on-chip op¬ 
tical isolators, unidirectional waveguides, and circula¬ 
tors [M]. Non-reciprocal bands predicted by topological 
band theory usually appears as surface modes attached 
to two-dimensional (2D) or three-dimensional (3D) bulk 
photonic systems. To make the device more compact, 
approaches based on symmetry breaking in waveguide 
structures usually suggest complex geometries such as 
helical structures [3IIM2]. 

Spectral reciprocity, cj(— k) = cj(k), can be protected 
by time reversal symmetry (T) and spatial symmetries 
such as inversion (V) [3J0], in addition to the local sym¬ 
metries in permitivity or permeability tensors (e T = e 
or /x T = pt). It is easy to understand that T symme¬ 
try can be broken by external static magnetic field [5], 
while V symmetry can be broken by using asymmetric 
structures such as chiral structures mmm or a sym¬ 
metrical structure under external magnetic field of spe¬ 
cific orientation m ■ However, spectral reciprocity can 
also be protected by a combination of symmetries such as 
spatial-temporal symmetries, which add more complexi¬ 
ties in the design of non-reciprocal waveguides. 

In this paper, we use compact non-chiral magnetized 
plasmonic waveguides consisting of only spherical parti¬ 
cles to demonstrate how spectral reciprocity can be pro¬ 
tected by a rotation-time-reversal ( 1ZT ) symmetry (i.e., 
time reverse followed by rotation of 180° about propaga¬ 
tion x-axis). In the 7 ZT symmetry broken case, we show 
that asymmetric dispersion relation can be supported. 
By coupling the hybridized bands with light lines, this 
simple system further supports one-way wave propaga¬ 
tion and energy transmission within a finite range of fre¬ 
quencies. 


one unit cell 



FIG. 1. (Color online) Geometry of a diatomic chain system. 
The chain contains two metallic nanoparticles in a unit cell, 
denoted by A and B. They can be elliptical or spherical, 
dependent on the values of semi-major axis a x and semi-minor 
axis ay, where a = A or B. In the case of spherical, a x — 
a^. A and B are separated with distance t y . External static 
magnetic field B is applied in z direction. Length of a unit 
cell is d. The chain breaks T, V, and 1ZT symmetry. 


II. MODEL AND METHODS 


We start by considering a magnetized diatomic sys¬ 
tem as shown in Fig. [l] The chain contains two types of 
metallic nanoparticles with different sizes, namely parti¬ 
cle A and B , with the same dielectric constant e(uj)/e o- 
The two nanoparticles formed “atoms” of a unit cell, and 
hence it is regarded as a dimer chain. As long as the 
nanoparticles are not too close together, the electromag¬ 
netic responses of the nanoparticles can be modeled by 
electric dipoles M- We denote the dipole moment of 
nanoparticle in the nth unit cell as p n;cr where a = A or 
B for type A and B particles, respectively. These dipole 
moments satisfy a set of self-consistent equations, known 
as the coupled dipole equations mm- 

^ ^ (^*-<7 ^nm^aa' ^Anmcrcr ') Pm;cr' — F n . a , (1) 

m,a' 

in which m runs from —N to N, a’ — A or B, 5 nrn is 
Kronecker delta function, and E®* 1 is the external driving 
field. We note that ol g is the quasi-static polarizability 
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FIG. 2. (Color online) Dispersion relations for an infinitely long and magnetized diatomic chain system without 1ZT symmetry 
(a) and with 1ZT symmetry (b). Sub-figure (i) is the chain geometry, while (ii) is corresponding dispersion relation obtained 
by density plot of 1/ |A|. (a)(ii) shows non-reciprocal bands (c o{k) ^ cj(— /c)), while only reciprocal bands (w(k) = cu(—k)) could 
be seen in (b)(ii). Bands related to z component are separated and not shown here. The blue solid lines are dispersions of 
free photon modes. Inner particle separation in (a) is t y = 0.75d, and in (b) is t x — 0.425^, where d! — 2d. Furthermore, 
(a)(ii) can be understood as hybridization of 4 bands, which are schematically drawn in dashed for guidance. Each dashed 
line represents different oscillation modes at k = 0, which is labeled by arrows at the middle. Ac o is the range at which only 
one-way propagation modes are allowed. 


with radiation correction of nanoparticle a and G nmcrcr / is 
interaction between dipoles p m;cr / and p n - a . Expressions 
are given in Appendix A [see Eqs. (A-l) and (A-4)]. 

We first consider the case without external driving 
field (E® x * = 0) and N —)> oo. Since the system 

is spatially periodic, by Bloch’s theorem we can write 
Pm;cr = Pk-,ae lkrnd , where k is wave vector. This simpli¬ 
fies Eq. 0 into a 6 by 6 matrix form (see Eq. ( | A-3| ) ): 


M k (u) 


Pk-,A 

Pk]B 


= o, 


where 


Mfc(a;) 


~ a/ O' 


GkAA GkAB 

i 

o 

p 

to ^ 


_ GkBA GkBB _ 


( 2 ) 


(3) 


In the above, Gkaa' = Go m(7(7 'e* fem(i , is the interac¬ 
tion between nanoparticles a and a' in k space. Instead 
of solving detM/ c (cj) = 0 to obtain the dispersion rela¬ 
tion psunnn], we may apply the eigen-response theory 
to evaluate the dispersion relation mm- Similar to 
the eigen-response theory, we plot 1/|A| as a function 
of k and cj, where |A| = |A(fc,u;)| is the smallest abso¬ 
lute value of the eigenvalue of the matrix M/ c (cj). This 
quantity gives huge value when there is resonance, and is 
plotted in Fig. (2|a). Note that dynamic dipolar Green’s 
function [lOj is used, and the infinite series in the inter¬ 
action up to \m\ = 120. 


III. FORMATION OF NONRECIPROCAL BANDS 


In Fig.^a), we showed the geometry and correspond¬ 
ing dispersion relation in (i) and (ii). This case consid¬ 
ers nanoparticles are spherical with radius a£ = 0.35d 
and a^ = ^0^5inner particle separation t y = 0.75d, 
and plasma wavelength A^ = c/(2 ttuj p ) = 10d, where uj p 
and c are plasma frequency and light speed in vacuum. 
Cyclotron frequency u) c = q\B\/m = 0.005co> p , in which 
B, m, and q are external static magnetic field, electron 
mass, and electron charge. As uo c oc |B|, it is treated 
as a variable to indicate the magnitude of B. Further¬ 
more, for simplicity, simple lossless Drude model is used 
m, and bands related to the £ component of pk - a , 
are not shown. It is because B is in z direction, so p^. a 
are coupled with p y k . a but not p z k . a , hence the bands can 
be separated. Note that the two solid blue lines are light 
lines, the dispersion of free photon modes. Region within 
light lines is light cone, modes in light cone are radiative 
and therefore not sustainable mmm- 

Figure [2ja) shows a case with non-reciprocal (asym¬ 
metric) bands. For comparison, using the same formal¬ 
ism, a case with reciprocal (symmetric) bands is shown 
in Fig. (2|b). This is a case where nanoparticle B and A 
are on the same axis, and horizontal inner separation is 
t x = 0.425d', where d! is length of unit cell in (b). We 
set d' = 2d, twice than that in (a), therefore the light 
cone in (b) is bigger than (a). Both (a) and (b) share 
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the same nanoparticles A and F>, cyclotron frequency cj c , 
and plasma wavelength (A p = 5 d! in this case). 

The non-reciprocal bands in Fig. |2ja) (ii) predicts only 
guided modes with positive group velocities are allowed 
within the range A uj. This results one-way propagation 
behavior, which can be utilized as an isolator. Note the 
operation frequency A uj is relatively broad, for example, 
about 50 times wider than the structure suggested in Ref. 
m- Range of operation frequency is about 1 x 10 a uj p in 
Ref. [TO] , while we have about 5 x lO -3 ^. The isolator 
also has a lower requirement on the external magnetic 
field, where cyclotron frequency uj c ~ 0.005a;p, is about 
10 times smaller than that used in Refs. uni na, and 
about 100 times smaller than that in Ref. [22] . 

We notice that the non-reciprocal bands in Fig. |2|a) 
are obtained by means of the simultaneous violation of 
certain symmetries, V, T, and 1ZT symmetries, whereas 
the only breaking of V and T still makes the system 
reciprocal shown in Fig.[2|b). We will discuss the relation 
between reciprocity and related symmetries below. 


(a) 


O O O O O ygv 

ooooo ® 

- External s 


External static 
magnetic field B 


OOOOO 


OOOOO ® B t 
L -^ 

= R 

Toooo 0..® B 
ooooo ® 


. Oo 

OO 



E T 

- Oo 

| R 

OO 

OO ■ ® B 

A 

u 

*oo 

Oo 

OO ® B 



FIG. 3. 1ZT acts on the diatomic chain systems shown in 
Fig. T is time reversal operation, which flips the direction 
of B when acting on the system. This can be understood 
as the motion of electrons are reversed EH- TZ is pi-rotation 
about x-axis. It exchanges the positions of A and B and flips 
the direction of B in case (a), only flips B in case (b). 1ZT is 
the operation T followed by 7 Z. For (a), as the final system is 
not identical to the initial one, the chain does not have 1ZT 
symmetry. For (b), the final system is identical to the initial 
one, the chain has 1ZT symmetry. 


from all nanoparticles in the infinite chain system [25]. 
So the T symmetry will be broken by external B only. 


A. Reciprocity protected by V and T symmetry 


B. Reciprocity protected by 1ZT symmetry 


Given p n;cr (£) = p fc;cr e z/cnd_2u;t is a solution of a system 
with frequency uj and wave vector k. Spatially inverted 
state will be V[p n](T (t)] = V {pk-a) e~ tknd ~ luJt . V turns 
(x,y,z) into (—x, — y, — z\ flips the direction of vector 
quantities, but does not modify uj. If the system has V 
symmetry, the inverted state will also be the solution of 
the system. The factor e - lknd ~ luJt Q f the new solution 
means that it is a solution with frequency uj but wave 
vector —k. This tells us that we are always able to find 
a solution with frequency uj but wave vector — k if the 
system has V symmetry, and thus the dispersion must 
be symmetric (uj(k) = uj{— k)). 

Time reversed state is obtained by taking complex con¬ 
jugate of the frequency component (see Ref. [23], also 
appendix E), 

T[PnAt)} = (Pk;ae iknd )* e~ iUt 

_ * —iknd—ibjt ( a\ 

— P/c;cr e • W 

This will be another solution if the system has T sym¬ 
metry. Similarly, it is a solution with frequency uj but 
wave vector — fc, as it has the factor e - lknd - luJt . Using 
arguments as those in V symmetry, we know the bands 
are symmetric as long as there is T symmetry. We see 
that, from Fig. [3] direction of external B will be flipped 
if T is operated on the system. Thus, the presence of 
external B breaks T symmetry, as the transformed one 
is not identical to the initial one. 

Note that lossy system will also break T symmetry. 
We assumed the material is lossless, and the effect of ra¬ 
diation loss is compensated if we count the contributions 


7 Z is pi-rotation operator, which rotates the system 
with 180° about x-axis. In the case of diatomic chain 
system shown in Fig. [l] the direction of external B is 
flipped, and positions of A and B are exchanged, see 
Fig. [3] 1ZT is the operation T followed by 1Z. When 1ZT 
acts on the solution, we have 

KT\p n .At)] = K{p% ia e- ilmd - iut ) 

= H(p* k;o .)e- iknd - i “ t , (5) 

here we used the property that 7 Z is not related to posi¬ 
tion x and time T. Again, if the system has 1ZT symme- 
try, Eq. (§ will be a solution with frequency uj but wave 
vector — k. This shows that the bands are symmetric 
about k. 

The matrix representation of 7 Z depends on system ge¬ 
ometry. For the chain system shown in Fig. [2ja) , 7 Z not 
only exchanges the positions of A and F>, but also rotates 
the vector quantities when acting on the states. Denot¬ 
ing the 3x3 rotation matrix which rotates a vector about 
x-axis with 180° by R = diag(l, —1, —1), then the 1ZT 
transformed state is: 


TZT [P n;a(t)\ 


' 0 R' 

’ Pfe;A ’ 

R 0 

. P t;B . 


iknd—iut 


(6a) 


For the chain system shown in Fig. |2](b), 1Z does not 
exchange positions of A and F>, so the 7 ZT transformed 
state is 


KT [Pn;<rW] 


R 0 

" Pt,A " 

0 R 

. Pfc;B . 


g— iknd—iujt 


(6b) 
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For the system of Fig. [2|a), the external B breaks 
the T symmetry, and simultaneously the non-identical 
nanoparticles A and B break the V and 7 ZT symmetry, 
see Fig. [3j The non-reciprocal band could then be ob¬ 
tained shown in Fig. [2ja) (ii). In contrast, when 1ZT 
acts on the system shown in Fig. [2^b) , B is flipped twice 
and remains unchanged. Meanwhile, particles A and B 
are on the x-axis, 7 Z would not modify their positions, 
and thus the transformed system is identical to the non- 
transformed one, which means it has 7 ZT symmetry, we 
get reciprocal dispersion relation in Fig. j2jb)(ii). 

Therefore, in order to obtain nonreciprocal bands in 
ID magnetized chain system, it is essential to break all 
related symmetries, including V, T and 7 ZT symmetries. 


IV. SYMMETRY OPERATORS ON DIATOMIC 
CHAIN SYSTEM 

A system is said to have 0 symmetry if it is invariant 
under transformation, 0 _1 M/ C (co’)0 = Mfc(o;). Here we 
show explicitly that the T and 7 ZT operator on the cou¬ 
pled dipole equation M/^cj). For simplicity, we consider 
the system shown in Fig. [lj employ quasi-static dipolar 
Geen’s function, use simple lossless Drude model, and 
neglect the radiation term in polarizability. 


A. Coupled dipole equation in quasi-static limit 


Quasistatic expressions are obtained by taking free 
space wave vector fco 0, so the polarizability, from 
Eq. (A-1), is 


Oi~ X = 


1 


eoV 


-uj 2 iujuj c 0 

-iujuj c —l d 2 0 I 


, (7) 


0 0 -cj 2 

in which we are assuming lossless model, with 7 = 0 . 
Symbols are defined under Eq. (A-l). Also, the quasi¬ 
static dipolar Green’s function, from Eq. (A-4), is 

1 

47re 0 


G ' k ™' - fc 0 7 o Gw - E c ^'( m ) 


Akrnd 


m^O 


r 3 

' 0 merer' 


, (8) 


where for the case shown in Fi g. [T| r elative position vec¬ 
tors romcrcr' are defined by Eq. (A-6), and 


2 0 0 

C AA (m) = C BB(rn) = I 0 -1 0 


/ 3 m 2 d 2 1 

m 2 d 2 +t 2 1 


C AB (m) = 


0 0-1 

3 mdty 
n 2 d 2 -\-t 2 


0 \ 




3 mdty 
m 2 d 2 -\-t 2 

0 


K 

m 2 d 2 -\-t 2 

0 


- 1 0 


- 1 / 


C BA(m) = C AB (-m) 


(9a) 

(9b) 

(9c) 


The above was obtained by putting Eq. (A-6) into (A-5) 


in appendix B. Some properties of G kaa , 
in appendix D. 


are discussed 


B. T and 1ZT transformation on the system 


When T operates on the system, it turns k into —k 
and takes complex conjugate (see also appendix E): 

T- 1 M fc (w)T = M_ fe (w)*. (10) 

In the quasi-static limit, Eq. ( flOj ) can be written as 

r’K- 1 - G' kaa ,)T = (ay 1 )* - G' w , (11) 


Therefore, T -1 M/ c (u;)T = Mfc(a;) only if (a^. -1 )* = 
a^ _1 , which means uj c = 0. That is, the system has 
T symmetry if external B = 0. 

For the 1ZT transformation, first we notice that R~ 1 = 
R, then (7^T) _1 M/ c (cj) 7^T = T7£Mfc(a;)7£T. In the 
quasi-static limit for the case shown in Fig. [lj we have 


TlZoc' 1 5 aa r'JZT 


= T 
= T 


J <7 <7 

0 R 

R 0 


a 


' - 1 0 
0 a B _1 


a 


! — 1 * 


a 


/ —l* 


T = 



'OR' 


R 0 


a 


T 


B 1 0 

0 ay 1 


(12a) 


and similarly, 


TTZG' kaa ,nT=G' k 


koo '' 


(12b) 


Eq. (12) implies that the 1ZT transformed system is not 


generally identical to the non-transformed one. So the 
chain in Fig. [l] has 1ZT symmetry only if a! A = a^, 
which is not true. 

The dispersion relation can be obtained by solving de¬ 
terminant equation det [M/ C (<y)] = 0. Neglecting z com¬ 
ponents, Appendix D shows that it is a polynomial with 
G'kAB.xy U P 2 n d order, and G' kAB xy is the only term 
which is odd in k in the determinant polynomial. Non¬ 
reciprocal dispersion only comes out when the polyno¬ 
mial is not an even function of &, which means coefficient 


of G[ 


kAB,xy 


, given by 


—2iujuj c 


1 

Va 


1 

Vb 


x (G' k 


G'u 


kAA,xx^kAB,yy 


GkAA,yyGkAB,xx) 


is non-zero. It is the case that both uj c ^ 0 and Va 7 ^ Vb, 
that is, external magnetic field B / 0 and the particles A 
and B are not identical. Therefore, T and 7 ZT symmetry 
should not be present in our case of Fig. 1. 
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V. ONE-WAY WAVE PROPAGATION AND 
ENERGY TRANSMISSION 

The non-reciprocal bands in Fig. [2ja)(ii) predicts only 
guided modes with positive group velocities are allowed 
within the range Acj, which gives one-way propagation 
behavior. We demonstrate the propagation behavior and 
energy transmission of the magnetized diatomic chain by 
considering it’s finite version in this section. 


A. One-way wave propagation 

For a finite magnetized diatomic chain system con¬ 
taining V unit cells, Eq. 0 can be written in a ma¬ 
trix equation form M(cj)p = E ext , where M(cj) is a 
( 12 V+ 6 ) x (12V + 6 ) square matrix vectorized from 
-M-nmcrcr' — &-<j (^) ^nm^aa' ^nmaa' • P &nd E are 
column vectors vectorized from p m;a / and E® x [., where 
each has (2V + 1) x 2 x 3 = 12V + 6 elements (2V + 1 
unit cells, 2 atoms per unit cell, and 3 spacial dimensions 
in our system). Since M(cj) is known, and E ext depends 
on our choices, so the excited dipole moments p can be 
found by evaluating the inverse: 

p = M(cj) - 1 E ext . (13) 

We study the finite version of magnetized diatomic chain 
in Fig.[ 2 ja), with V = 46 and plasma collision frequency 
7 = 0. Two types of driving polarizations are applied 
only to the site n = 0 , coherent in x direction, or coherent 
in y direction, see Fig. [4jb)(i) and (ii). In both cases, 
E e n f a = 0 for n ^ 0, while = (1,0,0) T for 

case (i), and E^ = Eg*g = ( 0 , 1 , 0 ) T for case (ii). 

The norm of excited dipole moments by three driving 
frequencies are shown in Fig. |4ja). Sub-figures ( 1 ) and 
( 2 ) correspond to coherent driving in x direction and y di¬ 
rection. From Fig. 0 a), we see the system supports two- 
way propagation at uj = 0.555cj p , as there are excitations 
throughout the chain; and supports one-way propagation 
at uj = 0.546(47,, as only spheres on the right are excited; 
and no supported modes at uj = 0.538cOp, as there is no 
excitation on both sides. It is reasonable, as slope of the 
dispersion relation is group velocity of the coupled plas- 
mon mode. From Fig. [2ja), we see at uj = 0.546^, only 
mode with positive k is allowed, therefore only mode with 
propagation is supported at this frequency. Note that 
reversing the propagation direction can be easily done by 
just flipping the direction of the external static magnetic 
field, this gives us a switchable optical isolator. 


B. One-way energy transmission 

Energy transmission is usually hard to define in plas- 
monic waveguides. Here we infer the transmission of the 
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FIG. 4. (Color online) (a) Dipole moments excited on a fi¬ 
nite magnetized diatomic chain system (93 cells). An exter¬ 
nal driving electric field is applied to the middle of the chain 
(n = 0) with driving frequency (i) u = 0.555 47 , (ii) 0.546^7 , 
(iii) 0.538u;p, are indicated by dotted lines in Fig. ga). The 
sub-figures in each frequency are corresponding to two types 
of driving polarization, coherent in x or y direction, labelled 
by (1) and (2). Note that at uj — 0.546c<7, one-way propa¬ 
gation occurred, as only modes with positive k are allowed, 
(b) External driving electric polarization, (i) coherent in x 
direction; (ii) coherent in y direction. 


diatomic chain shown in Fig. [l] by reading the dipole mo¬ 
ments Pn-a excited on nanoparticles. Since energy den¬ 
sity is proportional to square of electric field |Re(E)| , 
and since dipole moment satisfies p = aE, time aver¬ 
aged energy density of a particle a at cell n is thus pro¬ 
portional to |pn;cr | 2 = Pn;cr * Pn-cn which gives the sense 
of energy transmission. We consider the quantity defined 

by 


1 it-r** 

(|Pn| 2 ) = g E (lPm;A | 2 + |Pm;s| 2 )- 

m=n —4 

In the above we picked 4 cells near the nth for spatial 
average. One can improve by picking more cells, but 
in that case the chain has to be longer, or this is no 
longer a local quantity at cell n. This quantity is dis¬ 
played in Fig. |5j There are 4 driving polarization modes 
to excite the system, two in phase and two out of phase 
[26] , denoted by mode 1 to 4 as shown in the sub-figures 
of Fig. [5| It shows that there is no energy transmis¬ 
sion within band gaps, since {|P 35 1 2 ) = (|P- 35 | 2 ) = 0 at 
uj ~ 0.54co> p and uj ~ 0.61 uj p . The one-way property can 
be seen at uj = 0.546co> p in Fig. [5] (a)-(c), in which red 
curve is finite while blue curve is zero, implying there is 
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Energy transmission to the left or right ends (|p„ = _ 35 
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atomic chain can be understood as a result of four bands 
hybridization, colored in blue, yellow, red, and green in 
Fig. [2ja) (ii). We can see the transition of the four bands 
in Fig. [HJ 
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FIG. 5. (Color online) Energy transmissions to the left end 
(n s= —35) or the right end (n — 35). Energy density is 
proportional to the quantity (|p n | 2 ) (defined in text), and is 
plotted in the figure. Left (right) end quantity is denoted by 
the blue (red) curve and uses the left (right) upper horizontal 
frame ticks. The magnetized diatomic chain dispersion rela¬ 
tion is drawn as a background for reference. Sub-figures (a)- 
(d) are corresponding to different external driving polariza¬ 
tions, which are pictured in the sub-figures (mode 1-4). One¬ 
way propagation property is demonstrated at c o ~ 0.546^, 
at where red curve is finite but blue curve is zero. A little 
material damping 7 = O-OOfMceip is added to the nanoparticles 
here, which is to reduce extreme fluctuations so that excita¬ 
tions can be seen clearer. 


energy transmission to the right end but no transmission 
to the left end. 

From Fig. [5] we also see that different modes excite 
different frequency range. Polarization modes 1 to 4 cor- 

0.56co> p , and 


respond to uj > 0.55c uj < 0.58cu p , uj < 
uj > 0.62 uj p with respectively. This can be explained by 
using band hybridization model (discussed in next sec¬ 
tion), which is shown in Fig. [l a )(ii)- 


VI. HYBRIDIZATION OF BANDS IN DIATOMIC 
CHAIN 

There are four surface plasmon resonant modes for iso¬ 
lated dimer particles, two transverse modes (inphase and 
antiphase oscillations) and two longitudinal modes [27] . 
The four modes form dispersion bands when the dimer 
particles are duplicated, becoming a diatomic chain. Af¬ 
ter then, the dispersion relation of the magnetized di- 


FIG. 6 . (a) Dispersion relations of diatomic chains formed 
by different ellipsoids. Corresponding unit cell configurations 
are shown in (b). Moving from (i) to (vi), semi-minor axis a y 
are increased one by one, so nanoparticles are becoming more 
and more spherical. B = 0 , a x — 0.35d, and Vb/Va = 0.5 
in (i) to (vi). From (i) to (vi), the lower two bands move up¬ 
wards, cross the upper two bands, and hybridize each other. 
Note that (iv) is the case that A and B are spherical, (v) 
is different from (iv) by adding external static magnetic field 
(cj c = 0.005c o p ), resulting a non-reciprocal dispersion. In (vi), 
A and B are identical spheres with a x = of = 0.35d. Al¬ 
though B / 0, the dispersion is reciprocal as it is protected 
by V symmetry. 


For better understanding of the hybridization, here we 
consider the diatomic chain formed by ellipsoids. Spheres 
in Fig. [I] are replaced by ellipsoids with varying a x and 
a^, where t y = 0.75d. Fig. |6ja)(i)-(iv) show the dis¬ 
persion relations with increasing a y . In these 4 cases, 
a x = 0.35d, Vb/Va = 0.5, and B = 0. Corresponding 
unit cell structures are depicted in (b)(i) to (b)(iv). Fig. [6] 
(a)(i) shows four separated bands. The lower two bands 
move up and cross the upper two bands when the parti¬ 
cles are more spherical, as shown in (i)-(iv). If one further 
applies external static magnetic field to case (iv) such 
that uj c = 0.005^, we have case (v), which is the case in 
Fig.[2ja). Thus, in case (i)-(v), we see the deformation of 
bands, and so conclude that the dispersion in Fig. (6^iv) 
is formed by hybridization of the four bands. Case in 
Fig. [6jb)(vi) shows the case with nanoparticle B is iden¬ 
tical to nanoparticle A , that is a a x — a* = 0.35d, where 
a = A or B. This is obtained by replacing all particle 
B by A in Fig. [6jb)(v). Bands in Fig. |6jb)(vi) are sym¬ 
metric, as it has both 7 ZT and V symmetry. From these 
we know the exitance of B and non-identical spheres are 
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essential to achieve non-reciprocal bands. 

The hybridization model also explains the excitations 
in Fig. [5j Fig. [2ja)(ii) shows the upper two bands con¬ 
tain parts that are original from the red dashed band, 
therefore the excitation by polarization mode 1 is promi¬ 
nent at the upper two bands, as found in Fig. Up). This 
is because the oscillation mode of the red dashed band 
at k = 0 is the same as mode 1. Similarly, the lower 
two bands contain parts that are original from the green 
dashed band, and hence excitation by polarization mode 
2 is prominent at the lower two bands. 


where e = y 1 — a° 2 /a° 2 . As we consider nanoparti¬ 
cles A and B are of the same shape in this paper, so 
L a = L b- The last term with fcg accounts for the radia¬ 
tion correction [30] , which will be vanished in quasi-static 
approximation. 

With Drude model, the dielectric tensor is in the form 

mm 


e(u) _ ( £xx Exy 

— I ~ e xy e yy 

e ° \ 0 0 



(A- 2 ) 


CONCLUSIONS 

To conclude, we used a compact non-chiral magnetized 
plasmonic chain to demonstrate the crucial role of 7 ZT 
symmetry in the design of these subwavelength waveg¬ 
uides with non-reciprocal dispersion (cj(—k) 7 ^ u(k)). 
The hybridization among four plasmon modes and free 
photon modes give rise to a frequency range where only 
guided modes in one direction are allowed. While we are 
considering a weaker magnetic field, this operation fre¬ 
quency range is already much wider than that in Ref. 
m- Matrix representations of the symmetry operators 
were used to explain explicitly how the spectral reci¬ 
procity is protected by T, V, and 7 ZT symmetries. The 
results may provide a clear guidance on the design of 
one-way waveguides. 
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APPENDIX 


A. Polarizability and Drude model in external B 


The inverse polarizability of an ellipsoidal particle a 
(a° > dy = a°) is given by [10, 28, [29] 


OL 


-1 

<J 


1 

^V 



y|_j 

67re 0 


(A-l) 


where volume V a = 47 ra£a^a^/ 3 , I is identity 3x3 
matrix, fco is light wave vector in free space, and pa¬ 
rameter L a depends on particle’s shape, which is given 
by L a = diag (N x , N y , N z ). For spherical particles, 
N x = N y = N z = 1/3. For ellipsoidal particles, N x = 
(1 + e 2 )[ln i±f - 2e]/(2e 3 ) and N y = N z = (1 - N x )/2, 


in which e xx = e yy = 1 -w 2 (w + ry)/[w(w + ry) 2 -ww 2 )], 
e xy = -iw 2 w c /[w(w+* 7 ) 2 -ww 2 ], = l-w 2 /(w 2 +i 7 w). 

7 , and uj c are plasma collision frequency, plasmon 
frequency, and cyclotron frequency. In the case without 
material loss, A = 0, and e(u) becomes real. More terms 
have to be added to the static polarizability if one wants 
an accurate numerical evaluation on lossy materials, but 
it is not of interest in this paper. Further discussions can 
be found in Ch .8 in m- 


B. Coupled dipole equation in k space 


In Eq. (|T|) , defining !VI nmcrcr / cx (T S nrn S a(7 ' G nmaa 1 
, putting p m-a' = P k-a'^ lkrnd : and in the case E® 3 ^ = 0 , 
we have 

Xb ( J2 m M ~' eikmd ) Pfc;-' = 0- 

Multiplying both sides by e~ lknd : and notice M nmcrcr / 
depends onm-n only, we have M nmcrcr / = Mo ?m - n ,aa', 
and thus 

(E m M °.—^ ,e * fc(m n ^) P/e;a ' = 

As the sum m runs from — oo to oo, we have 

Writing it in matrix form, we have 


E 


° 

0 o. B 


^ 0 m 


Go rnAA Go mAB 

Akmd l 

Pfc;A 

Go mBA Gq mBB 

6 / 

_ P/c;£ _ 


0 . (A-3) 


This gives Eq. ©• 


C. Dynamic dipolar Green function for diatomic 
chain system 

The dipolar coupling between the particle a’ in the 
mth cell and the particle a in the nth cell depends on 














m - n only, i.e., G nrna(T > = G 0 , m - n ,aa', so we only 
show elements Go m(7(7 / here. Also, position vector of the 
particle a' in the rath cell is denoted by r mcr /. Rela¬ 
tive position vector is then ro m(JC r' = ro a — r mcr /, and 
the corresponding unit vector is therefore p a a'{m) = 
r 0 maa' /romaa' - Spatial components of the unit vector 
are denoted by p a a'(m)x, Paa'{m)yi and paa'{m)z- The 
dynamic coupling is well known and is given by uni 


G 


o ik 0 r 0 , 


Omaa' — 


47Te 0 
T Gq-q-' ijrri) 


/ \ ^0 
f[m)- 


^ 0 merer 

ik 0 


Omcrcr' 


Omcrcr' 


(A-4) 


for m 7 ^ 0 together with a' ^ cr, otherwise Goocrcr = 0 , as 
a particle is not interacting itself by generating electric 
field. In the above, eo is free space permittivity, and 
ko = uj/c is the light wave vector in free space. Matrixes 
in the above are 


/ (ra) = 

Paa'(m) 2 y 

~ Per cr' (jrri) x P cr cr' ( rn )y 


0 


Paa' { rn ) y Paa' { rn ) x 0 

Pcrcr' ( m ) x 

0 


and 

G era' {m) = 


^Pcrcr' (jTl)x 


1 


3 Paa' ( m) y p aa / (m) x 3/w {m) y - 1 

0 


3 P crcr' (m) x pa<r'(m) 
2 
y 

0 


y 


(A-5) 

For the system shown in Fig.[l]or[2ja)(i), relative position 
vectors are 


^*0 mAA — 1*0 mBB — ( md , 0, 0) , 
1*0 mAB = (-md , - ty , 0) T , 
and r 0mS A = {-md, t y , 0 ) T . 


(A- 6 ) 


For the system shown in Fig. §b)( i), relative position 
vectors are 


I*0mAA — 1*0 mBB — ( md , 0, 0) , 
romyis = (-md - t x , 0 , 0 ) T , 
and romBA = {—md, + t x , 0 , 0 ) T . 


(A-7) 


D. Properties of quasi-static Green function G' kaa / 
for system in Fig. |2|(a) 

There are some properties about the matrix G' kaa , re¬ 
garding to the system in Fig. [lj From Eq. (| 8 |, one can 


G'_ w = G', 


Also, from Eq. (9a) we know G' kAA = G' kBB , and they 
are both diagonal, so we have 


R- 1 G^R = Gi 44 = Gl 


kAA 


T kBB 


R 1 G' kBB H — G kBB — G kAA , 


(A-9a) 


where R = diag(l, — 1 , — 1 ) is the rotation matrix about 
x-axis with 180°. Furthermore, Eq. (9c) tells us that off 


diagonal elements in G' kAB and G' kBA are purely imagi¬ 
nary, and Eq. (9b) implies G kAB = G^ba , so we have 


R _1 Gi 


kAB 

R —lr-if t > c* r * 

— A — Vt I 


R = Gl 


kAB 

/ 

kBA 


= Gi 


kBA 

/ 

kAB • 


(A-9b) 


Elements in G' kaa , are functions of k. By expanding 
the terms and Eq. (A-9), we see that G' kAA , G' kAB xx and 
G'kAB yy are pnrely real and even in k , while G' kAB xy is 
purely imaginary and odd in k. Hence, we deduced that, 
neglecting 2 components, 


T kAA 


G'kAB ~ 


G'kBA — 


/ -( 

G'kAA,xx 0 

kBB ~ l 

0 G’kAA, 

c 

^ kAB,xx 

G'kAB,xy j 

rn 

^ kAB,xy 

G kAB,yy ) 

G'kAB,xx G kAB xy \ 

~G kABxy G kAByy J 


Therefore, the determinant 


\M k (u>)\ = 


OL’ 


L M-G 


-G'kBA 


kAA 


a 


GkAB 

/ — if, A n’ 


l B L ( CJ ) ^kAA 
is a polynomial with G' kABxy up to 2 nd order. 

E. The T operator 

Regarding to Eq. Q, noticing T _1 = T, we have 
^T _1 M nm Tp m = J2™ n mTp k e l 

m 

= Tj2 M n m vle~ 

m 

= Te~ iknd e iknd M nm p* k ( 

^—ik(m—n)d rt * p —iknd 

= TM-kPk e ~ iknd = M*_ fcPt e ftnii . 


Akmd 


—ikmd 


—ikmd 
u G 




In the above, we used the fact that M nm = Mo ?m _ n 
and defined M& = Mo m e^ md . The last line proves 


kaa' 


(A- 8 ) Eq. ( 10 ). 
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F. The 1ZT operator on diatomic chain system in 
Fig.gb) 

Double chain in Fig. |2^b) has 1ZT symmetry. Here we 
show that M/ c (cj) is commute with 7 ZT in quasi-static 
limit. By Eqs. ( |6b| ) and 


OL„ 


= T 


"R 0 " 


_ ° R_ 



a 


(RT) 

/ -l 
A 

0 


/ -1 


Oi 


Oil 


a 


= a'' -'S., 



"R 0 " 


°R 


T 


Also, by putting Eq. (A-7) into (A-5), we see C crcr /(m) = 
diag(2, —1, —1), and thus G' kaa , are diagonal. As a re¬ 
sult, R^Gfccrcr/R = GkcraA and we have 


C RT)~ l G 


kaa‘ 


(11T) = G, 


kcra' 


G. Comments on Pi-rotation ( 1Z ) symmetry and 
reflect ion-in-y ( V y ) symmetry 

As we assumed magnetic field is in ^-direction, we are 
actually considering the reciprocity of bands related to 
x and y components (abbreviated as xy bands), as 2 re¬ 
lated components are separated. The conclusion that 
“reciprocity is protected by 7 ZT symmetry” thus can 
be narrowed to “reciprocity is protected by 0T symme¬ 
try”, where 0 is the x-y plane projection of 7 Z. How¬ 
ever, reflection-in-^ (V y ), which changes (x,y,z) into 
(x, —y,z), has the same form under the projection. Ex¬ 
plicitly, matrix representations of the operators are, 

/1 0 0 \ / 1 0 0 \ 

R = 0 -1 0 and P y = 0 -1 0 , 

\0 0 -1 ) \0 0 1 / 

from which we see they are the same if we chop away 
z related components. As a result, 0 is corresponding 
to two spatial transformations, and the transformations 
have the same effects on xy bands. As the two transfor¬ 
mations are equivalent under the projection, so we may 
also arrive the conclusion that reciprocity is also pro¬ 
tected by T y T symmetry [32] . 
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